Electronic structure and effects of dynamical electron correlation in ferromagnetic 

bcc-Fe, fcc-Ni and antiferromagnetic NiO 
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We have constructed a LDA+DMFT method in the framework of the iterative perturbation theory 
(IPT) with full LDA Hamiltonian without mapping onto the effective Wannier orbitals. We then 
apply this LDA+DMFT method to ferromagnetic bcc-Fe and fcc-Ni as a test of transition metal, 
and to antiferromagnetic NiO as an example of transition metal oxide. In Fe and Ni, the width 
of occupied 3d bands is narrower than those in LDA and Ni 6 eV satellite appears. In NiO, the 
resultant electronic structure is of charge-transfer insulator type and the band gap is 4.3 eV. These 
results are in good agreement with the experimental XPS. The configuration mixing and dynamical 
correlation effects play a crucial role in these results. 
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I. INTRODUCTION 

Much attention has been focused on strongly corre- 
lated electron systems, since these systems show anoma- 
lous physical properties such as various spin, charge and 
orbital order, metal-insulator transition and so on. The 
local density approximation (LDA) based on the den- 
sity functional theory (DFT) is hardly applicable to these 
fruitful physical properties in strongly correlated electron 
systems. LDA+U method could not treat metallic phase 
near metal-insulator transition since it has been devel- 
oped in order to discuss magnetic insulators. 1 Fluctua- 
tion of charge or/and spin densities plays an important 
role in low-energy excitations near the metal-insulator 
transition. To discuss realistic materials of strongly cor- 
related electron systems, another method is needed for 
the "dynamical" electron correlation. 

The dynamical mean field theory (DMFT) has been 
developed and applied to model systems of strongly cor- 
related electrons, which led us to a unified picture of 
low- and high-energy excitations in anomalous metallic 
phase near metal-insulator transition. 2 DMFT is based 
on mapping of many electron systems in bulk onto single 
impurity atom embedded in effective medium, namely 
the single impurity problem. In this mapping proce- 
dure, the on-site dynamical correlation is included and 
the inter-atomic correlation is neglected. The combina- 
tion of DMFT with LDA, called LDA+DMFT, has been 
developed in order to discuss realistic systems. 3 ' 4 LDA 
or DFT meet their own new stage with the GW approx- 
imation, 5 which is based on the many body perturba- 
tion theory and the random phase approximation (RPA) . 
The combination of DMFT with the GW approximation 
(GW+DMFT) was proposed to include both on-site and 
inter-atomic Coulomb interaction. 6 

To solve the mapped single impurity problem within 
DMFT, one can use several computational schemes 
such as the quantum Monte Carlo method (QMC), 7-11 



the iterative perturbation theory (IPT), 12 15 the non- 
crossing approximation (NCA) 16 and the exact diagonal- 
ization (ED). 17 

Hirsch-Fye QMC 18 is widely used as a solver for the 
mapped single impurity problem within DMFT. How- 
ever, it is not applicable in the low temperature limit. 
Moreover, Hirsch-Fye QMC has serious difficulty in 
application of multi-orbital systems with spin-flip and 
pair-hopping terms of the exchange interactions, since 
one cannot apply the Hubbard-Stratonovich transforma- 
tion to these systems. 9 To go beyond such problems 
of Hirsch-Fye QMC, several QMC methods have been 
proposed. Projective QMC 10 has been developed to 
calculate in the low temperature limit and Continuous 
Time QMC (CTQMC) 11 to calculate multi-orbital sys- 
tems with spin-flip and pair-hopping terms. 

In spite of development of those QMC, the exact cal- 
culations such as QMC and ED can only be suitable 
for simple Hamiltonians with relatively small size be- 
cause of their computational costs. Several LDA+DMFT 
methods with QMC adopt the projected effective Hamil- 
tonian of Wannier-like functions with a few adopted 
bands for reducing computational costs. With the use 
of the effective Hamiltonian of Wannier-like functions in 
LDA+DMFT, only a few adopted bands are discussed 
with fixing the hybridization mixing and could not de- 
scribe a possible change of hybridization due to local 
Coulomb interaction. 

To carry out calculations in realistic materials with 
multiple-orbitals, approximate calculation schemes of 
IPT or NCA could be more suitable because of its ef- 
ficient CPU-time, though IPT is not applicable to cases 
of large Coulomb interactions and NCA cannot yield the 
Fermi liquid behavior at low energies and in low temper- 
ature limit. 19 

The IPT method was developed by Kajueter and 
Kotliar for non-degenerate orbital and the results show 
good agreements with the ED results. 12 However, it 
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was reported that the quasiparticle peak in the IPT re- 
sults tends to be sharper 2 and that the IPT results in 
Lai-^Sr^TiOa does not reprduce quasiparticle peak at 
1000K. 20 

We generalized IPT for multi-orbital bands and ap- 
plied this method to doubly degenerated c g bands and 
triply degenerated t 2g bands on simple cubic and body- 
centered cubic lattices. 13,14 We also verified that this 
DMFT with IPT scheme is fairly applicable to arbitrary 
electron occupation cases with different Coulomb interac- 
tion U. The spectrum shows electron affinity and ioniza- 
tion levels with different electron configurations of differ- 
ent occupation numbers and the system becomes insulat- 
ing state at an integer filling in sufficiently large U. 13 ' 14 
IPT for multi orbitals was also proposed by Laad, Craco 
and Muller-Hartmann and applied to several realistic ma- 
terials. 15 

In this paper, we apply this DMFT with IPT scheme to 
realistic materials, where we adopt full LDA Hamiltonian 
without reducing its size and IPT as a solver for the 
mapped single impurity problem. In the following, we 
notify this DMFT with IPT as "LDA+DMFT-IPT(1)" . 
The goal of the present paper is as follows: 

(i) To generalize LDA+DMFT-IPT(1) which is appli- 
cable to various realistic strongly correlated materials, 
both metallic and insulating, multi-atom (compound), 
spin-polarized and strongly hybridized cases between s, 
p and d-bands. 

(ii) To apply LDA+DMFT-IPT(1) to ferromagnetic 
bcc-Fe and fcc-Ni and to antiferromagnetic NiO. For fer- 
romagnetic bcc-Fe and fcc-Ni, we will discuss whether 
LDA+DMFT-IPT(1) reproduces the accurate width of 
occupied 3d bands and the observed satellite of spec- 
trum at 6 eV below the Fermi energy in Ni ("6 eV 
satellite"). Application of LDA+DMFT-IPT ( I ) to fer- 
romagnetic bcc-Fe and fcc-Ni is a test in transition met- 
als since previous LDA+DMFT 21 - 22 has been applied to 
those metals and reproduced those physical properties 
well. For antiferromagnetic NiO, we will discuss whether 
LDA+DMFT-IPT(I) reproduces the accurate band gap 
and the correct description of the charge-transfer insula- 
tor. These physical properties are not well reproduced in 
the previous LSDA calculation. 

(iii) To verify the applicability of LDA+DMFT- IPT(I) 
to various realistic strongly correlated materials. Ac- 
tually, the physical properties mentioned in (ii) are re- 
produced fairly well by LDA+DMFT-IPT(1). We will 
discuss the origin of drastic changes of those physical 
properties and show the validity of LDA+DMFT-IPT ( I ) 
by comparing with other theoretical scheme including 
other LDA+DMFT method 21-23 and experiments. We 
will discuss the origin of nickel 6 eV satellite, which 
has been concluded as hole-hole scattering process, 24 by 
using the spectrum of Ligand Field Theory (LFT) ob- 
tained in the procedure of LDA+DMFT-IPT(1). In bcc- 
Fe and fcc-Ni case, the narrowing of occupied 3d bands in 
LDA+DMFT-IPT(1) occurs due to the on-site dynam- 
ical electron correlation. In NiO case, the appearance 



of charge-transfer insulator and accurate band-gap in 
LDA+DMFT-IPT(1) is due to the change of hybridiza- 
tion between Ni-3d and 0-2p bands caused by local 
Coulomb interaction of nickel 3d bands. 

In Sec. II, the general formulation of LDA+DMFT- 
IPT^) will be given. In Sec. Ill, we test LDA+DMFT- 
IPT^) in systems of ferromagnetic bcc-Fe and fcc-Ni. 
In the energy spectra, one can see nickel 6 eV satellite 
caused by hole-hole scattering process and the narrow- 
ing of occupied 3d bands caused by on-site dynamical 
electron correlation. Section IV will be devoted to cal- 
culation of the electronic structures in antiferromagnetic 
NiO by LDA+DMFT-IPT ( I ) . Local Coulomb interac- 
tion of nickel 3d bands enhances the hybridization be- 
tween Ni-3d and 0-2p bands and system becomes the 
charge-transfer insulator with accurate band gap. In 
both sections, we will present the energy spectra and k- 
resolved spectrum and LFT spectra of single isolated ion. 
Section V is the summary. We discuss the applicability 
of LDA+DMFT- IPT(l) to various realistic strongly cor- 
related materials by comparing with other LDA+DMFT 
method and experiments. 



II. FORMULATION 

A. Hamiltonian and Coulomb interactions 

We proceed with the Hubbard-type Hamiltonian as fol- 
lows : 

H = #LDA + Hmt ' W 

^LDA = ^2 (^mm'cr ^ _ Ah% c 6jj> S mm > )c] m(y Cj' m' a , (2) 
j(ma)j'{m') 

Hint = ^ ^ ' ^TOim 2 :TO3TO4 C jm 1 CT 1 C jm2<T2 C J m 4CT4 C jni3<T3 ' 

j(mcr) 

(3) 

The Hamiltonian H is divided into two parts: the 
unperturbed part (Hlda) an< ^ the interaction part 
{Him). ^w° A) in Eq. (2) is the full LDA Hamiltonian 
constructed by the Tight-Binding LMTO (TB-LMTO) 
method, 25 without projecting onto any kind of effective 
local orbitals. The use of original full LMTO Hamilto- 
nian in LDA+DMFT-IPT(1) enables us to describe the 
change of hybridization caused by local Coulomb inter- 
action. Ahj c in Eq. (2) is the double-counting term in- 
cluded in LDA Hamiltonian as averaged Coulomb and 
exchange interactions, because we need to subtract A/i^ c 

from h^J^^^ in order to construct the unperturbed part 
of the Hamiltonian Hf^ A . Note that an index j in the 
summation of Eqs. (2) and (3) runs over atomic sites, 
{m} orbital indices and {a} spins. 

U mi m 2 :m 3 mi hi Eq. (3) is the on-site electron-electron 
interaction matrix and presented by using the Slater in- 
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tegrals F k as 



^ a k {mim 2 : m 3 m 4 )F k ,(4:) 

a<k<2l 

k 



Oi(mim 2 : m 3 m 4 ) 



= 2k+ \ ^ ( lm l\ Y kq\ lm 2l 
q— — k 

x (/m 3 |y fe * 9 |Zm 4 ), 



(5) 



where {|Zm)} are the basis set of complex spherical har- 
monics. Note that Eqs. (4) and (5) are the general ex- 
pression and invariant under any rotational operation. 1 
For d-orbitals we only need F°,F 2 and F 4 . The av- 
eraged Coulomb and exchange parameters U and J for 
3d-orbitals are defined in terms of the Slater integrals 
F°, F 2 and F 4 . Here we regard these equations as a def- 
inition of the Slater integrals F°, F 2 and F 4 in terms of 
U and J: 



U = F , 

J = j~ A (F 2 + F 4 ), 
with a ratio of F 2 and F 4 for 3d orbitals as 26 



F 4 /F 2 



0.625. 



(6) 
(7) 

(8) 



Therefore, once we obtain experimentally or assume the 
values U and J of 3d transition metal elements, the values 
of the Slater integrals F°, F 2 and F 4 can be evaluated by 
using Eqs. (6), (7) and (8) and then each matrix element 
of the on-site electron-electron interaction U, 



can be evaluated by using Eq. (4) and (5). 

As discussed above, we treat the interaction Hamil- 
tonian Eq. (3) for full 3d-orbitals on transition met- 
als within the framework of DMFT. Thus, the present 
Hamiltonian in Eq. (3) does work for the cases with par- 
tially filled e g and t2 S -orbitals of ferromagnetic bcc-Fe, 
fcc-Ni and antifcrromagnctic NiO. On the other hand, 
the case of triply degenerated t2 9 -orbitals and doubly de- 
generated Hubbard model were treated within the frame- 
work of DMFT in Rcfs. 15 and 27, respectively. 



B. Dynamical mean field theory 

The matrices of the lattice Green's function and the 
local Green's function are presented on the basis of the 
non-orthogonal local base set as 

^(fe,^)]" 1 

= [(i Wn + /x)O(fc) - {H^ A (k) + S(iw n )}], (9) 



G(ioj n ) = J dk G(k,iu n ), 



(10) 



where H^ A (k), 0(k) and fj, are the Hamiltonian matrix 
of Eq. (2), the overlap matrix, both in the fc-space, and 
the chemical potential, respectively. Here we neglect the 



fc-dependence of the self-energy E within the framework 
of DMFT. The chemical potential \i is determined to sat- 
isfy Luttinger's theorem. 28 All the matrices in Eqs. (9) 
and (10) have suffices {jma,j'm'a'}, where j, m and a 
stand for atomic sites in a unit cell, orbitals and spins, 
respectively. 

We assume that the Green's function in the effective 
medium G°(iu) n ) may be expressed as the self consistent 
equation of DMFT: ' 



G°(iu> n ) 1 =G(iuj n ) 1 + E(iw n ) + p, - /i, 



(11) 



where jl is the chemical potential in the effective 
medium. 



C. Iterative Perturbation Theory 

In this paper, IPT method is adopted as a solver for 
the mapped single impurity problem, generalized by Fuji- 
wara et al. for multi-orbital bands on arbitrary electron 
occupation with different Coulomb U. 13,14 In the IPT 
scheme, the self-energy is determined with the interpo- 
lation scheme between the high frequency limit and the 
strong interaction limit (the atomic limit) by using the 
second order self-energy. 

The second order self-energy is calculated as 



£< 2 ) (t)= 



i<W) 



x U mi m 3 :mm 2 G miCy (T)G t l l3Cr , (T)G° n2Cr / (— t) (12) 



vj( 2 )(^„)=/ dre l ^ T ^(T) . 



(13) 



By using the second order self-energy £( 2 )(^„), the self- 
energy is then expressed in a matrix form as 

S(iw„) = S HF + {iuj n ) [1 - B(iu n )XM (iu n )] -i , (14) 

where E HF is the Hartree-Fock energy. The matrices A 
and B(iui n ) are determined by requiring the self-energy 
Yi{iu) n ) to be exact in the high-frequency limit (iu> n ~ * 
oo), and in the atomic limit (U — > oo) 13 as 

A= lim {E(iw„)-E ifF }E( 2 )(iw„)- 1 , (15) 



lLO n — >OQ 



B(ioj n ) - E^ 2 )^)" 1 - {E at (^„) 



}- 1 A,(16) 



where E at (iw„) and E a *( 2 )(iw n ) in Eq. (16) are the self- 
energy and the second order self-energy in the atomic 
limit. These values are calculated exactly within the 
framework of the ligand field theory (LFT), which is dis- 
cussed in Section II D. 

The self-energy in Eq. (14) also includes spin-flip and 
pair-hopping terms of the exchange interactions since the 
second order self-energy in Eq. (12) includes those terms. 
In addition, the IPT can be achieved within efficient CPU 
time since the self-energy in Eq. (14) is obtained by only 
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matrix calculation and hence the IPT is applicable to re- 
alistic materials with larger size by using full-LDA Hamil- 
tonian. Form the discussion in this subsection, present 
IPT is much appropriate impurity solver and present IPT 
appropriately treats on-site dynamical correlation effects 
including both hole-hole and electron-electron scattering 
as well as electron-hole scattering. 

In the following Sections (Sees. Ill and IV), we assume 
the cubic symmetry around transition metal ions and, 
therefore, the self-energy can be assumed as 



i a >{iu) n ) — 'Ej mcr (iu! n )5 mm '5jj'S 



(17) 



For a cubic symmetry, the local Green's function becomes 
diagonal after the fc-integration. 



D. Calculation in the atomic limit based on the 
ligand field theory 

The Hamiltonian mapped onto an isolated atom is de- 
fined as 



Hatom 

Tjdc 
n atom 



Tjdc I TT 

n atom ~r n int , 

^ ^ €mm' o-cl na -C m f a , 
(ma) 

, (LDA) r 
n mm'a ~ ^n dc Omm' ■ 



(18) 
(19) 

(20) 



The first term H^ om in Eq. (18) is an unperturbed one- 
electron part. The Coulomb interaction Hint is the same 
as Eq. (3). /^nrnv i s tne on -site element of LDA Hamil- 
tonian and therefore ft^*' includes the information of 
bulk matrix with electron transfer, crystal field and or- 
bital hybridization. Ah dc is the double counting term 
same as in Eq. (2). Then one electron energy t-mm'a in 
Eq. (20) includes the crystal field splitting. Once we ob- 
tain e mm 'a and the Slater integrals, many-electron eigen- 
states of Eq. (18) should be obtained within the frame- 
work of the ligand field theory (LFT). 29 Full configura- 
tion interaction (CI) calculation is carried out by using 
the basis of all the multiplet of d electrons, the number 
of which is 2 10 = 1024. Using Full-CI calculation yields 
an exact solution for local self-energy of Eq. (14) in the 
atomic limit. 

The Green's function, the self-energy and the second- 
order self-energy in the atomic limit are then defined as 



Gmm'A^n) 



^mm'a( iuJ ) 



{vlCmal^iMn'M 



Z iu n + (E n - (j,Nr,) - {E v - ijlN v ) 

x { e -/3(£„-/iiV„) +e -l3(E v -^N v )yp^ 
ICO - (emm'a - - [G°* Mlw, . (22) 



mm2:mim 3 



m\ni2ms<y' 



x U, 



mim 3 :mm 2 



frri2cr' (1 fm\o fm^a' ) ~^ fmlO■fm30^ , 



(23) 



where (77), E ni N v , f ma are the many-electron eigenstates 
in the Hilbert space of 2 10 multiplets, the energy eigen- 
value, the total electron number of the eigenstate \r}) and 
the Fermi distribution function obtained from full-CI cal- 
culation within the framework of LFT. 

One can assign, with a help of the atomic spectrum 
in LFT, the origin of peaks in LDA+DMFT spectra, 
e.g. the initial and final multiplets of the corresponding 
transitions. Generally, peak structures in LFT spectra 
would suggest the existence of strong scattering channels 
in LDA+DMFT spectra at around these particular ener- 
gies. 

From discussion in Sections II C and II D, the IPT is 
a appropriate approximation method and is worth be- 
ing extended to the application of realistic materials in 
following four reasons: (1) Efficient computational cost 
(2) Applicability of realistic materials with larger size 
and more complicated hybridization by using full-LDA 
Hamiltonian [This is followed by (1)] (3) Care of spin-flip 
and pair-hopping terms of the exchange interactions in 
multi-orbital systems (4) Direct comparison of LFT spec- 
tra within the IPT method with the LDA+DMFT spec- 
tra and assignment of the origin of peaks in LDA+DMFT 
spectra 



E. Energy spectrum and fc-resolved spectrum 

The energy spectrum and the fc-resolved spectrum 
may be more physically understandable if we present the 
Green's function in the orthogonal base, subject to the 
sum rule of the energy integration of the local Green's 
function. The lattice Green's function and the local 
Green's function in the orthogonal base are defined as 

[S?(fc, w)]-^ 0~i (fc)G- 1 (fc, Lo)O-^(k) (24) 

-0-i(k){H^ A (k) + ^)}0-i(k)], (25) 
Sf(w) = i J dk <g{k,uj). (26) 

The energy spectrum will be shown with the imaginary 
part of fc-resolved spectrum will be also pre- 

sented with an imaginary part of Sf(fc,w), which may 
be compared with the angle-resolved photoemission spec- 
tra (ARPES). 

F. Computational Details 

The 8x8x8 fc-points are used in the fc-integration 
within the whole Brillouin zone. The fc-integration 
is carried out by using a generalized tetrahedron 
method. 13 The Pade approximation is adopted for ana- 
lytic continuation of the Green's function from the Mat- 
subara frequencies iu> n to the real w-axis. We adopt 
2 11 = 2048 Matsubara frequencies. 
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Here, all the calculation were carried out at T = 
1000 K to reduce the number of adopted Matsubara fre- 
quency. To consider temperature dependence of mag- 
netic order transition, one should include entropy cal- 
culation as well as total energy calculation. 30 Here, we 
do not include entropy calculation within LDA+DMFT 
scheme. Thus temperature T = 1000K ~ O.leV causes 
the spectrum broadening by at most 0.1 eV and that do 
not change the magnetic order. 

We first carry out LDA calculation to obtain input 
LDA Hamiltonian and then carry out DMFT calcula- 
tion. We adopt the converged DMFT results as output. 
We do not feed back the density by DMFT to the LDA 
calculation. 

In NiO case, on-site self-energy of nickel 3d bands 
on each sublattice is only included and inter-atom self- 
energy between nickel atoms on different sublattices are 
neglected. Thus, our calculation is "single-site" DMFT, 
but not cluster-DMFT. 31 



III. FERROMAGNETIC BCC-FE AND FCC-NI 

LSDA calculation has achieved a great success on over- 
all understanding of the ground state electronic struc- 
tures of crystalline transition metals, including bcc-Fe 
and fcc-Ni. However, LDA overestimates the width of oc- 
cupied 3d bands for both ferromagnetic bcc-Fe and fcc-Ni 
and could not produce the observed satellite of spectrum 
at 6 eV below the Fermi energy in Ni ( "6 eV satellite" ) . 
On the other hand, LDA+DMFT calculation 21 ' 22 gives 
the width of occupied 3d bands for ferromagnetic bcc-Fe 
and fcc-Ni and the Ni 6 eV satellite, in good agreement 
with XPS experimental results. 32,33 

The aim of the calculation for ferromagnetic bcc-Fe 
and fcc-Ni in this section is to see whether LDA+DMFT- 
IPT(l) scheme can reproduce a reasonable width of oc- 
cupied 3d bands and Ni 6 eV satellite. In fact, we will 
show that LDA+DMFT-IPT(1) presents the width of oc- 
cupied 3d bands for bcc-Fe and fcc-Ni and the 6 eV satel- 
lite in fcc-Ni to be in good agreement with both experi- 
mental ARPES 32 - 33 and also with previous LDA+DMFT 
results. 21,22 



A. Hamiltonian and U, J values 

Lattice constants, atomic sphere radii in the TB- 
LMTO method and the averaged values of Coulomb and 
exchange interactions, U and J, for ferromagnetic bcc-Fe 
and fcc-Ni to construct the Hamiltonian are summarized 
in Table I. The values of U and J are obtained with 
the constraint-LDA calculation. 21 The basis set of the 
Hamiltonian consists of full valence bands of 4s, 4p and 
3d orbitals, totally nine orbitals. 



B. Energy Spectrum 

Figure 1 shows the energy spectrum — Im^(o>) ob- 
tained by LDA and LDA+DMFT- IPT ( 1 ) with the 
atomic spectrum by LFT. The occupied 3d-band becomes 
narrower in LDA+DMFT- IPT(l) than in LDA, narrower 
by 0.3 eV and 1.0 eV for bcc-Fe and fcc-Ni, respectively. 
Particularly, the spectra near the Fermi energy, in the 
region —1.5 eV < u < 0.0 eV for both bcc-Fe and fcc-Ni 
in LDA+DMFT- IPT(l), becomes narrower than that in 
the low energy region. These results are attributed to 
strong rcnormalization of the quasiparticle caused by the 
on-site dynamical electron correlation within the frame- 
work of DMFT. The precise width of occupied 3d bands 
is discussed in Sec. IIIC. 

In addition, the satellite appears in LDA+DMFT- 
IPT(l) spectrum of fcc-Ni at 6 eV below the Fermi 
energy, in good agreement with experimental XPS re- 
sults. 33,36 Moreover, this satellite structure is strongly 
spin-dependent, much enhanced in the spectrum of the 
majority spin. These effects mainly come from the mul- 
tiplet scattering of d 8 — > d 7 , which is assigned by spin- 
dependent peak of LFT spectrum of d 8 — ► d 7 at around 
—5 eV in Fig. l-(2c). The satellite is due to hole-hole 
scattering process, 24 and can appear only when one treats 
whole on-site dynamical correlation effects. It should be 
noted that 6 eV satellite is not observed in GW approxi- 
mation (GWA) 37 since GWA includes the dynamical cor- 
relation only within RPA (electron-hole excitations) and 
thus not the electron-electron and hole-hole scattering 
processes. 

Let us focus on the atomic spectrum obtained by LFT 
in Fig. 1. The atomic spectra are separately shown in 
Figs, l-(lc) (2c) for ionization spectra and Figs, l-(ld) 
(2d) for affinity spectra. The initial state is a mixture of 
d 6 , d 7 in Fe and d 8 , d 9 in Ni, since the total occupation 
numbers are non-integers, nf^ t = 6.829 and nf^ t = 8.799 
from the results of bulk systems. The excited state is 
then a mixture of d 5 , d 6 , d 7 , d 8 in Fc and d 7 , d 8 , d 9 , d 10 
in Ni. Small components appear in ionization spectra 
above Ep and in affinity spectra below Ep. Moreover, 
characteristic atomic spectra at around —0.5 eV in Fe and 
Ni originate from both ionization and affinity process; 
d 6 -> d 7 and d 7 -> d 6 in Fe in Figs. l-(lc)(ld) and 
d 8 -> d 9 and d 9 -» d 8 in Ni in Figs. l-(2c)(2d). These 
effects comes from the fact that higher energy occupied 
states are not fully occupied. 

Characteristic structure in both atomic ionization and 
affinity spectra originate from various multiplet scatter- 
ing; d 7 — » d 6 and d 6 — > d 5 at around -4.0 eV and 
-6.0 eV and d 7 -> d 8 and d 6 d 7 at around 1.5 eV 
in Fe of Figs. l-(lc)(ld) and d 9 -> d 8 and d 8 -> d 7 at 
around -5.5 eV in Ni of Figs. l-(2c)(2d). This is due 
to the multiplet scattering caused by the exchange in- 
teraction J. The primary splitting of multiplets comes 
from the Coulomb interaction U, for example the tran- 
sition spectrum of d 8 — > d 7 should locate in the lower 
energy region than that of d 9 — > d 8 of Ni spectra. Then 
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FIG. 1: Left (la)-(ld): Energy spectrum in ferromagnetic bcc-Fe. Right (2a)-(2d): Energy spectrum in ferromagnetic fcc-Ni. 
(la) (2a): LDA. (lb)(2b): LDA+DMFT-IPT(1). (lc)(2c): Atomic ionization spectrum by the LFT. (ld)(2d): Atomic affinity 
spectrum by the LFT. The energy zeroth is set at the Fermi level (Ef = 0), and temperature is set to be 1000 K (T = 1000 K). 
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TABLE I: Lattice constants a, atomic spheres radii so, the Coulomb and exchange interactions U and J for ferromagnetic 
bcc-Fe, fcc-Ni and antiferromagnetic NiO. 





lattice constant 


a (A) atomic sphere radius so (A) U (eV) 


J (eV) 


fcc-Ni 
bcc-Fe 
NiO 


o.oZoo 

2.8708 
4.1948 


1.3768 3.0^ 
1.4128 2.0 21 
1.2318 (Ni), 1.0699 (O), 0.8882 (ES) 7.0 34 


0.9 Z1 
0.9 21 
0.9 35 


TABLE II: 

that Wd.occ 


The magnetic moment 
is not shown in 22. 


/z S pi n and the width of occupied 3d bands VKd.occ for ferromagnetic bcc-Fe and fcc-Ni. Note 






LSDA GW 37 LDA+DMFT LDA+DMFT-IPT(1) 


Expt. 


A t spin(pB) 


bcc-Fe 
fcc-Ni 


2.17 2.31 2.28™ 2.16 
0.47 0.55 0.57 22 0.47 


2.13 38 
0.57 38 


W d , occ (eV) 


bcc-Fe 
fcc-Ni 


3.7 3.4 3.4 
4.5 3.3 3.5 


3.3 32 
3.2 33 




TABLE III: The band 


gap -Egap (eV) and the magnetic moment ^t sp in (ps) f° r antiferromagnetic NiO. 






LSDA LSDA+ 


U 35 GW 41 ' 42 QPscGW 43 LDA+DMFT 23 LDA+DMFT-IPT(1) 


Expt. 


S ga p(eV) 

Mspin(Ms) 


0.2 3.7 
1.00 1.59 


0.2 4.8 4.3 4.3 
1.00 1.72 1.85 1.00 


4.3 39 
1.64 44 



bcc-Fe fcc-Ni 




ft ft 

FIG. 2: The fc-resolved spectrum — InrSf (k, ui) by LDA+DMFT-IPT(1) (shaded regions) with the energy bands (dashed lines) by 
LDA. Left: Ferromagnetic bcc-Fe: (a) Majority spin, (b) Minority spin. The high-symmetry fc-points are F(0, 0,0), N(i, |,0), 
P(tt, |, |), H(0, 0, 1). Right: Ferromagnetic fcc-Ni: (a) Majority spin, (b) Minority spin. The high-symmetry fc-points are 
L(i, |, |), r(0,0,0), X(0, 1,0), W(|, 1,0), K(0, |, |). The energy zeroth is set at the Fermi energy E F = 0, and temperature 
is set to be 1000 K (T = 1000 K). 
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additional multiplet splitting is caused by J. In LFT 
spectra in —6.0 eV < u> < 1.5 eV, the energetic order of 
the spectral positions does not follow the above first sim- 
ple rule and this fact implies that the scattering process 
by J changes the multiplet spectra drastically. We con- 
clude that the structures in LDA+DMFT-IPT ( 1 ) spectra 
is broadened and smoothed due to the multiplet scatter- 
ing of d 7 -> d 6 and d 6 -> d 5 at -6.0 eV < u < -2.0 eV 
and that of d 7 — ► d 8 and d 6 — > d 7 at around 1.5 eV 
in bcc-Fe and that of d 9 — > d 8 and d 8 — » d 7 at around 
-6.0 cV < u < -2.0 eV in fcc-Ni. 



C. fc-resolved spectrum and magnetic moment 

Figure 2 shows the fc-resolved spectrum — Im^(fc,u;) 
by LDA+DMFT-IPT(1) with the LDA energy bands. 
LDA results overestimate the width of occupied 3d va- 
lence bands, which is defined to be the energy differ- 
ence between the Fermi energy and the energy eigenvalue 
at P-point (Fe) or L-point (Ni), in comparison with ex- 
perimental results. The width of occupied 3d valence 
bands and the magnetic moment in LDA+DMFT-IPT(1) 
are shown in Table I in comparison with those by LDA 
and the experiments. 32,33 ' 38 The valence band width of 
LDA+DMFT-IPT(1) is in reasonable agreement with ex- 
periments. 

One can observe flat branches in bcc-Fe at around 
— 1.5 eV of majority spin and that at around —2 cV of 
minority spin. These flat bands correspond to the local 
multiplet excitations of d 7 — > d 6 . One can also observe 
flat branches in fcc-Ni at —6 eV of both majority and 
minority spins. These flat bands at —6 eV are due to 
the local multiplet excitations of d 8 — ► d 7 , where the 
intensity of the fc-resolved spectrum of majority spin is 
larger than that of minority spin. This intensity differ- 
ence causes the strong spin dependence of — 6 eV satellite 
in the energy spectrum in Fig. l-(2b). 

The fc-resolved spectrum in bcc-Fe is more diffusive 
than that in fcc-Ni at around —3.0 eV ~ 0.0 eV in spite 
of smaller U for bcc-Fe than fcc-Ni. This comes from that 
the scattering process by J compared with U is stronger 
in bcc-Fe than in fcc-Ni since J/U in bcc-Fe is larger than 
that in fcc-Ni. 

The magnetic moment /i sp m of both bcc-Fe and fcc- 
Ni is almost the same as LDA result. /i sp i n of bcc- 
Fe in LDA+DMFT-IPT ( 1 ) result is also in good agree- 
ment with experiment, while that of fcc-Ni is slightly 
smaller than experimental result. This is an artifact due 
to insufficient number of fc-points to integrate the lat- 
tice Green's function G(k,ui) to obtain the local Green's 
function G(oj) by using generalized tetrahedron method. 
The adopted value of the total number of the fc-points in 
LDA+DMFT-IPT(1) is much smaller than that in TB- 
LMTO method and the discrepancy of ^ sp ; n for fcc-Ni 
may be improved by increasing the total number of the 
fc-points in LDA+DMFT-IPT(1). 



D. Comparison with previous LDA+DMFT results 

Here, we compare the results for ferromagnetic bcc-Fe 
and fcc-Ni obtained by LDA+DMFT- IPT(l) with other 
previous LDA+DMFT 21,22 . Lichtenstein et al. 21 has 
used the LDA+DMFT with QMC as an impurity solver. 
Minar et al. 22 has used the KKR+DMFT with pcrturba- 
tive SPTF (spin-polarized T-matrix+FLEX) as an impu- 
rity solver. The energy spectra in LDA+DMFT-IPT(1) 
show a good agreement with those LDA+DMFT calcula- 
tions. The existence of Ni 6 eV satellite is also very simi- 
lar to those LDA+DMFT calculations, presumably much 
better coincident with the position of the observed spec- 
tra. The magnetic moment of bcc-Fe is in good agree- 
ment with those LDA+DMFT results. Slightly smaller 
value of the magnetic moment of fcc-Ni in LDA+DMFT- 
IPT^) than previous LDA+DMFT is not due to the use 
of perturbative IPT approach as an impurity solver but 
due to insufficient number of fc-points mentioned above. 
Thus, we can conclude that LDA+DMFT-IPT(1) repro- 
duces reasonable results for ferromagnetic bcc-Fe and 
fcc-Ni and that LDA+DMFT- IPT(l) is applicable to re- 
alistic metallic materials in strongly correlated electron 
systems as well as other previous LDA+DMFT meth- 
ods. 21,22 



IV. ANTIFERROMAGNETIC NIO 

NiO is a type-II antiferromagnetic insulator with Neel 
temperature 523 K. The experimentally observed band 
gap is 4.3 eV. 39 Resonance photocmission experiments 40 
show that electronic structure of NiO should be of the 
charge-transfer type and the low energy satellite mainly 
consists of nickel 3d bands. 

Various theoretical methods have been applied to NiO. 
The band gap -E gap and the magnetic moment ^ sp i n are 
summarized in Table III. LSDA calculation shows that 
NiO is Mott-Hubbard type insulator with a small band 
gap of 0.2 eV for antiferromagnetic phase, 45 in which 
oxygen 2p bands are located at lower energy region than 
the occupied nickel 3d bands. 

The LSDA+U method was applied and the resultant 
band gap is almost good agreement with experiment and 
a system becomes charge-transfer insulator. 35 However, 
bonding states of nickel e g bands are observed at around 
8 eV below the Fermi energy and this causes less compo- 
nents of the occupied main peak and too much compo- 
nents of the occupied satellite peak, compared with ex- 
perimental XPS spectrum. These problems mainly come 
from the static potential correction with orbital depen- 
dence in LSDA+U and this implies that one should in- 
clude dynamical correlation effects. 

GW approximation was applied and the resultant band 
gap still remains very small (0.2cV). 41,42 Quasiparticle 
self-consistent GW (QPscGW) approximation 43 was also 
applied and a system becomes charge-transfer insulator. 
However, the resultant band gap is overestimated (4.8eV) 
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and the bonding states of nickel e g bands are observed 
at around 6 eV below the Fermi energy. These fact im- 
plies that the dynamical electron correlation should play 
a more crucial role in NiO than treated by RPA in GW 
approximation. 

LDA+DMFT was applied to paramagnetic NiO 23 ' 46 
and the resultant band gap is in good agreement with 
experimentally observed XPS result. 39 and the system 
becomes charge-transfer insulator. However, the top of 
valence bands mainly consists of nickel 3d bands, while 
the cluster-model CI calculation 34,47 and oxygen x-ray 
absorption of Li^Nii-^O 48 show that the top of va- 
lence bands is mainly based on oxygen 2p bands. This 
difference comes from that these LDA+DMFT calcula- 
tions use projected effective Hamiltonian constructed by 
Wannicr-likc functions with fixing the hybridization mix- 
ing. Moreover, no calculation has been carried out for 
antiferromagnetic NiO by LDA+DMFT. 

The target of LDA+DMFT calculation is to get elec- 
tronic structures of antiferromagnetic NiO with (i) the 
accurate band gap, (ii) the correct description of the 
charge-transfer type insulator and (iii) the occupied main 
peak corresponding to oxygen 2p states and the occupied 
satellite peak corresponding to nickel 3d states. 



A. Hamiltonian, the Coulomb and exchange 
interactions U and J 

The structure of NiO is not of the dense packing and we 
put empty atom spheres in vacant region of the lattice in 
LMTO formalism. The lattice constants, atomic sphere 
radii of each atom and the averaged values of Coulomb 
and exchange interactions, U and J, are summarized in 
Table I. The total number of atom spheres in an anti- 
ferromagnetic unit cell is eight; two Ni, two O and four 
empty atoms (ES). We use the full LDA Hamiltonian 
of the LMTO formalism and adopted muffin-tin orbitals 
arc 4s, 4p and 3d in Ni and 2s, 2p in O and Is, lp in 
ES and, thus, the total number of basis is 42. U refers 
to the experimental value of U = 7.0 cV 34 and J to the 
constraint-LDA calculation of J = 0.9. cV 35 The differ- 
ence of [/-values in fcc-Ni and NiO is due to the difference 
of screening mechanism in metals and insulators. 

The value of U obtained by constraint LDA is U = 
8.0 eV, 35 which is slightly larger than experimental value. 
Within the framework of constraint LDA, all the screen- 
ing channels are switch off. When we evaluate suit- 
able value of U for 3d-orbitals required in the frame- 
work of LDA+DMFT, only the screening channels asso- 
ciated with 3d electrons should be switch off and those 
associated with 4s and 4p electrons should remain. In 
this sense, a suitable value of U is slightly smaller than 
that for constraint LDA. Thus, the adopted value of 
U = 7.0 eV in the present paper is reasonable within 
the framework of LDA+DMFT. 



B. Energy Spectrum 

Figure 3 shows the energy spectrum — Im^(w) of anti- 
ferromagnetic NiO by LDA+DMFT-IPT(1) with exper- 
imental XPS spectrum. 39 Experimental XPS spectrum 
mainly consists of three parts: a main peak at 4 eV, a 
main peak at —3.5 eV < u> < —0.5 eV and a satellite 
peak at -10.0 cV < uo < -6.0 eV as in Fig. 3-(4). With 
the cluster model CI calculation, 34 these three structures 
are assigned to d 9 , d 8 L and d 7 final state, respectively, 
where L is a ligand hole created in oxygen 2p orbitals. 

The energy spectrum by LDA+DMFT-IPT(1) repro- 
duces these three structures fairly well, the positions of 
two main peaks are in good agreement with experimen- 
tally observed XPS result, but the position of satellite 
peak at —10.0 eV < ui < —6.0 eV shifts slightly upward 
in comparison with experimental results. 39 Appreciable 
component of oxygen 2p bands appears just at the top of 
the valence bands. The main peak of the occupied states 
at —1.5 eV and the satellite at —6.0 eV <~ —7.5 eV origi- 
nate from t 2 g orbitals with minority spin of Ni and c g +t 2g 
orbitals with majority spin of Ni, respectively. The main 
peak of conduction bands comes from the Ni-e ff orbitals 
with minority spin. The spectrum of LDA+DMFT- 
IPT(l) in Fig. 3 shows that the electronic structure is 
of the charge-transfer insulator type, while the electronic 
structure in LDA is of Mott-Hubbard type. In fact, the 
hybridization mixing between oxygen 2p bands and nickel 
3d bands is much enhanced in LDA+DMFT-IPT(1) in 
comparison with that in LDA, though the Coulomb ma- 
trix elements between Ni-3d and 0-2p bands and among 
0-2p bands are not included. 

On-site Coulomb interaction between nickel e g bands 
makes the occupied nickel e g bands shift to the lower 
energy side and the unoccupied nickel e g bands to the 
higher energy side. Due to the shift of the occupied 
nickel e g bands to lower energy side, the hybridization 
between Ni-e g and Ni-t2 S increases and the occupied Ni- 
t 2g bands shift to the lower energy side. Since the oc- 
cupied Ni-e g and t 2g bands shift to the lower energy 
side, the hybridization between Ni-3d and 0-2p bands 
increases and whole bands are broadened in the region 
—8.5 eV < uj < —0.5 eV. These effects enhance the 
hybridization mixing and change the character of NiO 
to be of the charge-transfer insulator type. The energy 
spectrum and the electronic structure by LDA+DMFT- 
IPT(l) is in good agreement with XPS experiments. 

Atomic spectrum for single nickel ion obtained by LFT 
is shown in Fig. 3-(3). Initial state is 3 A2 9 (d 8 ) of the 
electron configuration ti^t^e 2 ^. The excited states 

of atomic ionization spectra al and a2 are 4 Ti g (d 7 ), bl 
and b2 are 2 E s (d 7 ) and cl and c2 are 2 Ti ff (d 7 ), respec- 
tively. This assignment is perfectly consistent with pre- 
vious LFT calculation. 34 

In the atomic calculation of LFT, a single nickel ion 
is considered and inter- atomic electron transfer is not al- 
lowed. Therefore, present atomic calculation does not 
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FIG. 3: Energy spectra — Im^(u;) of antiferromagnetic NiO. (1) Partial — Im5f(u) by LDA. (2) Partial — ImSf(t<;) by 
LDA+DMFT-IPT(1). (3) Partial -ImSf(tj) in the atomic spectra of Ni by LFT. (4) Total -ImSf(w) by LDA+DMFT-IPT(1) 
and XPS. The energy zeroth is set at the Fermi energy (Ef = 0), and temperature is set to be 1000 K (T=1000 K). Char- 
acteristic atomic spectra for e s (t2 S ) orbitals is labeled by white (black) letter in (3). Note that the index of nickel "majority 
(minority) spin" denotes up (down) spin component on up spin Ni site and down (up) spin component on down spin Ni site. 
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include explicitly the d L final state. The initial state 
in the LFT spectra is fixed in d 8 since the total elec- 
tron number of Ni ion is eight in NiO. Therefore, the 
unoccupied and occupied spectra in Fig. 3- (3) originate 
from the transition of d 8 — > d 9 and d 8 — > d 7 , respec- 
tively. Though the d 8 L final state is not included in 
LFT, of Fig. 3-(3), the occupied main peak in the spec- 
trum of LDA+DMFT-IPT(1) is in good agreement with 
experimental XPS result. In fact, the initial states in 
LDA+DMFT-IPT(1) is a mixture of d 8 and d 9 L, since 
the electron occupation numbers of nickel 3d bands and 
oxygen 2p bands in LDA+DMFT-IPT(1) are 8.17 and 
5.14, respectively. Therefore, we can assign, in the spec- 
trum of LDA+DMFT-IPT(1), the unoccupied main peak 
and occupied satellite peak to be the d 9 , d 7 final state, re- 
spectively. The occupied main peak is also assigned to 
be the d 8 L final state. 

The splitting between d 8 L and d 7 final state configu- 
rations is small and, in Fig. 3-(4), the occupied satellite 
peak in LDA+DMFT-IPT(1) appear at higher energy 
region (—6 ~ —7 eV) than that in experimental XPS re- 
sult (—8 ~ —9 eV). This may be due to the fact that 
we do not include the d 8 L final state in LFT calculation. 
The cluster model Cl-calculation instead of LFT calcula- 
tion would give more precise position of occupied satellite 
peak. 

Unoccupied e g bands of nickel majority spin locates at 
slightly higher energy region than that of minority spin 
and oxygen 2p bands in Fig. 3-(l). It would not be the 
case, if we include the inter-atomic component of the self- 
energy or cluster Cl-calculation. 

The band gap -E gap and the magnetic moment fi sp i n are 
summarized in Table III, The band gap in LDA+DMFT- 
IPT(l) is in fairly good agreement with experimental re- 
sult. On the other hand, the calculated magnetic mo- 
ment is almost the same as the calculated result by LDA 
and much smaller than the observed one. The magnetic 
moment of antiferromagnetic NiO comes from the elec- 
tron occupation of nickel c g bands since nickel t2 S bands 
are almost fully occupied. The present result of smaller 
value of the magnetic moment in LDA+DMFT-IPT(1) 
than the experimental result comes from that unoccu- 
pied e g bands of nickel majority spin has some intensity. 
In the present calculation, we neglect the inter-atomic 
and inter-spin components of the self-energy and the dis- 
crepancy may be attributed to this approximation, which 
should be left as a future study. 



C. fc-resolved Spectrum 

The fe-resolved spectrum — ImS^(fc,a;) by 
LDA+DMFT-IPT(1) with both the energy bands 
by LDA (dashed lines) and the angle-resolved pho- 
toemission spectrum (dots) is shown in Fig. 4. LDA 
band at 0.0 eV < u> < 1.5 eV is shifted to the 
LDA+DMFT-IPT(1) band at 3.0 eV < u < 4.0 eV. 
These spectra correspond to unoccupied nickel e g bands 



in Fig. 3-(l). These unoccupied nickel 3d bands in 
LDA+DMFT-IPT(1) may be more localized than that 
in LDA. 

In LDA band structure, the occupied bands mainly 
consists of two parts: (i) —2.5 eV < u> < 0.0 eV 
bands structure mainly based on Ni-3d bands and (ii) 
—8.5 eV < u < —3.5 eV bands structure mainly based 
on 0-2p bands. In LDA+DMFT-IPT(1) results, (i) and 
(ii) have been shifted to lower and higher energy side, re- 
spectively. This shift of these results causes the increase 
of hybridization between Ni-3d and 0-2p bands. Espe- 
cially, strong hybridization occurs at L and T points. 

We observe a broad and flat diffusive structure at 
-5 eV in LDA+DMFT-IPT(1) spectrum, which corre- 
sponds to occupied satellite peak of Ni-t2 9 bands and 
also Ni-t2 fl state in LFT in Figs. 3-(2)(3). Therefore, we 
conclude that Ni-t2 9 bands may be more localized due to 
the strong scattering by U. 

One can see flat branches at — 5.5eV < u) < 
—5.0 eVand at —6.5 eV < u < —6.0 eV in fc-resolved 
spectrum of LDA+DMFT-IPT(1). This corresponds to 
occupied satellite peak of Ni-e g bands in Fig. 3-(l). This 
flat band implies that the occupied satellite peak of Ni-e g 
bands is not due to the band structure but to the mul- 
tiplet scattering, particularly d 8 — ► d 7 assigned by the 
atomic calculation of LFT. 

The fc-resolved spectrum at —5.5 eV < u < —2.5 eV of 
LDA+DMFT-IPT(1) is in good agreement with experi- 
mental ARPES result 49 along T to X point, which is due 
to the inclusion of dynamical correlation within DMFT, 
which causes strong hybridization between Ni-3d and O- 
2p bands. 

Finally, we should mention that much enhanced inten- 
sity profile can be seen in -ImG(fe,w) but — Im^(fc,w) 
becomes more diffusive due to overlap integrals. 



D. Comparison with other calculations and 
experiments 

Here, we compare the results for antiferromagnetic 
NiO obtained by LDA+DMFT-IPT ( 1 ) with other pre- 
vious calculations. LDA+DMFT-IPT(1) reproduces 
4.3 eV band gap and the characteristics of the charge- 
transfer type insulator. Moreover, the results obtained 
by LDA+DMFT-IPT(1) shows no Ni-c 3 bonding states 
below the Fermi energy, which is observed in LSDA+U 35 
and QPscGW results. 43 LDA+DMFT-IPT ( 1 ) also repro- 
duces the top of valence bands with mainly based on 
oxygen 2p bands, while application of the LDA+DMFT 
method using Hamiltonian constructed by Wannicr-like 
function to paramagnetic NiO 23,46 shows that the top 
of valence bands is mainly based on nickel e g bands. 
Those effects obtained by LDA+DMFT-IPT(1) are in 
good agreement with experimental XPS results. Those 
drastic change in LDA+DMFT-IPT ( 1 ) is due to enhance- 
ment of hybridization between Ni-3d and 0-2p bands 
caused by on-site Coulomb interaction of nickel 3d bands 
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within DMFT scheme. 

The assignment of final states of the main and satel- 
lite occupied peaks and the unoccupied main obtained 
in LDA+DMFT spectrum is consistent with that in the 
cluster-model CI calculation. 34 This comes from the use 
of LFT spectrum within the IPT method and hence the 
assignment of the origin of peaks in LDA+DMFT spec- 
tra by using LFT spectra is a great advantage of the IPT 
method, compared with other impurity solvers. 

However we still have several problems, slightly higher 
energy positions of the occupied satellite peak and the 
unoccupied main peak and smaller value of the mag- 
netic moment. From discussion in the previous two 
subsections, these are mainly due to neglect of the 
inter-atomic and inter-spin components of the self-energy 
within DMFT scheme. Particularly, dynamical correla- 
tion of nickel 3d bands between two different sublattices 
and that of nickel 3d and oxygen 2p bands are not in- 
cluded since DMFT neglects all the inter-atomic dynam- 
ical correlation. 

To improve those problems in antiferromagnctic NiO, 
one may adopt cluster DMFT 31 , which is based on map- 
ping of many electron systems in bulk onto single clus- 
ter impurity problem with including the intra-cluster dy- 
namical Coulomb interaction and neglecting the inter- 
cluster Coulomb interaction. In LDA+DMFT-IPT(1), 
the extension of LDA+DMFT to LDA+cluster DMFT 
is consistent with the extension of LFT of single nickel 
isolated ion to the cluster model Cl-calculation of NiC>6 
cluster similar to Ref. 34 and 47. The use of the cluster 
model Cl-calculation will reproduce better positioning of 
occupied satellite peak and direct assignment of occupied 
main peak with d 8 L final state. 

From above discussion, we can conclude that 
LDA+DMFT-IPT(1) reproduces reasonable results for 
antiferromagnetic NiO within the DMFT scheme, and 
that LDA+DMFT-IPT(1) is applicable to realistic com- 
pound cases and insulating cases as well as LDA+DMFT 
with nonperturbative impurity solver. 

V. CONCLUSION 

In summary, we have proposed LDA+DMFT-IPT(1) 
where we use the full LDA Hamiltonian and IPT as an 



impurity solver. We then applied LDA+DMFT-IPT(1) 
to ferromagnetic bcc-Fe, fcc-Ni and antiferromagnetic 
NiO. 

For bcc-Fe and fcc-Ni case, on-site dynamical correla- 
tion effect causes the narrowing of the width of occupied 
3d bands. The multiplet scattering effect of d 8 — > d 7 
assigned by the LFT spectrum causes the spin depen- 
dent 6 eV satellite for fcc-Ni. The energy spectra and 
magnetic moment in LDA+DMFT-IPT(1) show a good 
agreement with previous LDA+DMFT calculations 21,22 . 

For antiferromagnetic NiO case, on-site Coulomb in- 
teraction of nickel 3d bands enhances the hybridization 
between Ni-3d and 0-2p bands and this causes the band 
gap of 4.3 eV, the charge-transfer type insulator, and 
the top of valence bands with mainly oxygen 2p bands. 
Those effects are all in good agreement with experiments. 
The drastic change of hybridization between nickel 3d 
and oxygen 2p bands is due to the use of full-LDA 
Hamiltonian and on-site dynamical electron correlation 
of nickel site within DMFT. The successful application of 
LDA+DMFT-IPT(1) to antiferromagnetic NiO implies 
that LDA+DMFT-IPT(1) is fairly applicable to larger 
compound cases with more complicated hybridization be- 
tween more than two atoms. 

In addition, the assignment of characteristic peaks of 
NiO in LDA+DMFT spectrum is consistent with that 
in the cluster- model CI calculation. 34 Thus, the assign- 
ment of the origin of peaks by using LFT spectra is a 
great advantage of the IPT method in understanding the 
existence of strong scattering channels in LDA+DMFT 
spectra. 

Several remained problems in antiferromagnetic NiO, 
positions of the occupied satellite and the unoccupied 
main peak, would be solved by including the inter-atomic 
dynamical correlation. 

Thus, we can conclude that LDA+DMFT- IPT(l) is 
applicable to various realistic materials, such as both 
metallic and insulating cases, multi-atom (compound) 
cases, spin-polarized cases and strongly hybridized cases 
between s, p and d-bands. 
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